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Abstract. The mass and anisotropy of an elliptical galaxy can be si- 
\& ' multaneously determined from velocity dispersion and line profile shape 

measurements. We describe the principles, techniques, and limitations of 
this approach, and the results obtained sofar. We briefly discuss how best 
to combine these stellar-dynamical results with X-ray measurements and 
Q\ ' gravitational lensing analyses. 

in ' 

1. Introduction 

OO ' 

ON ■ To place elliptical galaxies in the context of galaxy formation theories, it is vital 

^ e— | to understand their mass distributions and dynamical structure. Yet despite 

Qh' nearly two decades of work we are only beginning to answer important questions 

q ! such as: 

5— i ' > What are the density distribution, the flattening, the extent of elliptical galaxy 

c/2 . halos? Are these similar as for the halos of spiral galaxies, for the same stellar 

mass? 

[> To what extent has the luminous matter in ellipticals dissipatively segregated 
in the dark matter potential? Is there a "conspiracy" between luminous and 
dark matter to produce a flat circular orbit rotation curve, like in spiral galax- 
ies? 

> How anisotropic, radially or tangentially, are the velocity ellipsoids? 
[> How do the derived halo properties and dynamical structure compare with 
those obtained from numerical simulations of the merger of two spiral galaxies? 
Can the results be used to constrain the role of the baryons in hierarchical 
merging? 

Progress in this field has been beset by well-known difficulties: There is no 
single ideal tracer like a spiral galaxy's rotation curve, and the interpretation of 
velocity dispersion data is ambigous, suffering from a fundamental degeneracy 
between anisotropy and mass. However, recent advances on several fronts give 
reason for optimism: 

o From absorption line profile data it is now possible to determine both the 
mass distribution and anisotropy out to, at present, < 2 — 3 half-light radii 
R e . How this can be done, and what the results are to date, is the subject of 
this article. The method works, in principle, for all nearby ellipticals. 
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o The improved spatial and spectral resolution of ROSAT and ASCA has al- 
lowed accurate mass determinations to large radii for several massive ellipticals 
(Mushotzky et al. 1994, Kim & Fabbiano 1995). Resolution limits the analysis 
to R ?t R e . The number of X-ray bright ellipticals analysed in this way is still 
relatively small. 

> Galaxy-galaxy (strong) lensing gives accurate surface mass measurements at 
a specific radius defined by the lensed images (e.g., Kochanek & Keeton 1997). 
Galaxy-galaxy weak lensing (Brainerd, Blandford & Smail 1996) will in future 
allow a statistical determination of outer surface mass profiles for galaxies at 
typical redshift of a few tenths. 
Both X-ray and gravitational lensing mass measurements will be described else- 
where in these proceedings. In this article we explain the stellar dynamical 
principles of how line profile data can break the traditional degeneracy (Section 
2). We show that the analysis of real data does indeed work (Section 3), and 
then describe the results obtained sofar (Section 4). It is important to note that 
combining results from several approaches will be very useful; thus we briefly 
comment on how the new dynamical determinations fit in with the results from 
other techniques (Section 5). 

2. Mass and anisotropy from line profile data: Principles 

Illustrative line-of-sight velocity profiles (vps) for isotropic, radially, and tan- 
gentially anisotropic distribution functions in fixed potentials are given, e.g, in 
Dejonghe (1987), Merritt (1987) and Gerhard (1991) for spherical systems, and 
in Dehnen & Gerhard (1993, 1994), Evans (1993), and Qian et al. (1995) for 
axisymmetric systems. The signature of radial orbits is a broad central pro- 
file, and a peaked profile with extended wings beyond the half-mass radius, 
where the profile is dominated by orbits near tangent point, the point along the 
line-of-sight closest to the galaxy center. Tangential orbits tend to generate 
narrow profiles in the center and broad, sometimes double-peaked profiles at 
large radii. In axisymmetric potentials, rotation introduces a noticeable asym- 
metry with respect to the vp center, and the minor and major axis anisotropies 
can be decoupled. 

While the anisotropy of an elliptical galaxy is thus clearly constrained by 
vp shape measurements, it can not be directly inferred because the vps also 
depend on the gravitational potential (Gerhard 1993). An updated illustration 
for spherical Jaffe models of various anisotropy structure is given in Fig. The 
vp shapes are conveniently represented by Gauss-Hermite moments (Gerhard 
1993, van der Marel & Franx 1993). In the spherical case the vp is symmetric, 
and is parametrized in Fig. p] by the fitted profile width and lowest order even 
coefficient /14 (van der Marel & Franx 1993). Fig. |l] shows a general trend: As 
the mass of the model at large r is increased at constant anisotropy, both the 
projected dispersion and /14 increase. Increasing (3 at constant potential, on the 
other hand, lowers a and increases /14. This suggests that by modelling a and 
/14 both mass M(r) and anisotropy (3(r) can in principle be found. 

Locally, this trend may be reversed, as shown by the two Osipkov-Merritt- 
models in Fig. |l|; in this case a large number of high-energy radial orbits all 
turning around near radius r a lead to flat-topped vps in a small radial range 
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Figure 1. Anisotropy parameter (3, projected velocity dispersion a !it , and VP- 
parameter h 4 for representative Jaffe models in self-consistent (left) and halo poten- 
tial (right). CT f it and h 4 are obtained by simultaneously fitting a Gaussian and a 4 th - 
order Gauss-Hermite function to the VP. The models shown are the isotropic model 
(solid) , radially and tangentially anisotropic models constructed with the method of 
Gerhard (1991) (dashed lines), and two radially anisotropic Osipkov (1979)-Merritt 
(1985) models (dotted lines). Note that while /3 is a function of three-dimensional 
radius r and ant, h 4 are observed quantities depending on projected radius R, there 
is a close correspondence between features in these profiles. 



near r a . Although the properties of these models are extreme and they have not 
been very useful sofar for modelling observed vps, this nonetheless suggests that 
a global non-parametric analysis of the data is required. 



3. Analysis of observed line profile data 

In recent studies, analyzing simulated and real line-profile data has proceeded 
in two distinct steps (Merritt 1993, Rix et al. 1997, Gerhard et al. 1997): First, 
the gravitational potential <1> is held fixed, and in this current potential a best- 
estimate df is determined from the data. Then, this process is repeated for 
different potentials, and confidence intervals are found for the parameters enter- 
ing the functional form of 

In the first step, the df is written as a sum of elements (which can be values 
in grid cells, orbits, or basis functions). These elements are projected into ob- 
servational space, and their weights are adjusted to fit the data. The equation 
relating the df to the projected df (over line-of-sight velocities and sky positions) 
is linear, so this is relatively straightforward. Also, at least in the spherical case 
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we know that the solution is unique for ideal data (Dejonghe & Merritt 1992). 
The same is likely to be true for edge-on axisymmetric systems, as long as only 
one orbit family is involved. Because real data provide only a sparse and noisy 
subset of the projected df, only gross features of the phase-space df can be 
extracted, and significant smoothing is necessary. This is done by the concept of 
regularisation (see Merritt 1993). For each set of data an optimal regularisation 
parameter A must be determined. 

Fig. Q shows the application of such a technique to Monte Carlo data, which 
is obtained from a model df in a known potential, but using the sampling and 
radial extent of real data for the EO galaxy NGC 6703 (Gerhard et al. 1997). 
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Figure 2. (Left): Model analysis of pseudo-data generated from a radially 
anisotropic model DF. The sampling and radial extent are the same as for the real 
data of NGC 6703 (Fig. n), but assuming error bars Acr = 3kms _1 and A/i 4 = 0.0T. 
The full curves show the true <r, hi, and /3 profiles of the underlying DF. The 
dashed lines show the projected a and ft 4 profiles derived from the pseudo-data by 
our method for optimal A. The dotted curves show the uncertainties that remain in 
the deprojected quantity f3 even with such small error bars. 

Figure 3. (Right): Top: The cumulative distribution of the normalized xl+hi, 
for dynamical models recovered from 100 random Gaussian data sets derived from 
the same DF, with observational errors as in Fig. El (dashed line) and Fig. M (full). 
Bottom: Cumulative distribution of residuals between true and recovered distribu- 
tion function, evaluated on a grid extending to three times the radius of the last 
data point, for the same 100 datasets. Also shown in the top panel are the cumula- 
tive Y^ +h4 -distributions for data drawn from radially anisotropic mpdels in the two 
potentials corresponding to two of the extreme solid lines in Fig. |4[ with the same 
error bars as in Fig. |6| (dot-dashed lines) , and for a self-consistent model with more 
complicated anisotropy structure (dotted line). 
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The assumed sizes of the error bars in this example are 3kms for the velocity 
dispersion a, and 0.02 for the vp shape parameter /14. The elements used in 
this particular technique are various isotropic and tangentially anisotropic basis 
functions constructed along the lines described by Gerhard (1991). Fig. [2] shows 
that the true df can be recovered well from such data, and Fig. ^ shows quan- 
titatively the cumulative distribution of x 2 P er data point for random Monte 
Carlo realizations. In this example, the algorithm will match the data with 
Xa+h4 < 1-28 with 95% confidence, and recover an underlying smooth df to 
within an rms of 12%. 

The second step, that of constraining the potential, is harder because $ enters 
non-linearly into the equation for the projected df. Moreover, experiments sug- 
gest that the potential is in practice determined only nearly uniquely (Gerhard 
et al. 1997). In other words, there is likely a finite region in $-space within 
which for each $ there exists a positive df that matches even exquisite kine- 
matic data. It is not practical at present to optimize the gravitational potential 
non-parametrically. Thus in work sofar a functional form with one or two pa- 
rameters has been assumed, and confidence limits were obtained for the values 
of these parameters. See Section 4 for an example. Then one must, of course, 
assess whether the results are dependent on the chosen parametrization. 

4. Results sofar: The EO galaxy NGC 6703, and others 

In this Section we discuss the results obtained sofar. We will concentrate on the 
case of NGC 6703 which we have analyzed (Gerhard et al. 1997), but mention 
briefly the results obtained for NGC 2434 by Rix et al. (1997). Both galaxies 
are EO, so presumably nearly spherical, which simplifies the analysis. 

The kinematic data for NGC 6703 were obtained during three observing runs 
at the Calar Alto 3.5-m telescope with a Boiler & Chivens longslit twin spec- 
trograph, during a total of 18.5 hrs of observations. A major axis spectrum and 
two spectra perpendicular to the major axis and shifted by 24" and 36" from the 
center were taken. The instrumental resolution was 85kms . The data were 
reduced as described by Bender (1990) and Bender, Saglia & Gerhard (1994); 
extensive Monte Carlo simulations were used to test the effects of continuum 
fitting and instrumental resolution, and to determine the error bars. The data 
for a and /14 are shown in Fig. ^ below. 

Dynamical models were fit to these data along the lines described in the pre- 
vious section. For the luminous component we used a Jafle model which fits the 
surface brightness data to < 0.1 mag and the cumulative mass profile to < 2%. 
Models in which this luminous component has constant M/L and contains all 
the mass in NGC 6703, fail to fit the data by a large margin in x^+m (Fig- !)• 
Either the velocity dispersion profile is fit approximately, but then the line pro- 
files cannot be reproduced, or the /14 profile is fitted well, but then a(R) is poorly 
matched. The non-parametric algorithm finds a bad compromise which is ruled 
out at high statistical significance. 

Models with dark matter were investigated using a two-parameter family of 
dark halo potentials and rotation curves: 

M*0 = H Hr 2 + r 2 c ), v c {r) = *>0^=|- (1) 
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Figure 4. (Left): Rotation curves for a sequence of gravitational potentials 
(stars plus dark halo) used in the analysis of NGC 6703. The full lines show rota- 
tion curves that are consistent with the NGC 6703 kinematic data inside the 95% 
confidence boundary for A= 10~ 5 (open symbols in Fig. H). The other line styles 
show rotation curves inconsistent with the data; among these is the constant-M/L 
model with no dark halo (short-dashed). 

Figure 5. (Right): Quality with which the kinematics of NGC 6703 can be 
fitted in different potentials. The figure shows the average x 2 P er a ~ an d hi data 
point, of the best-estimate distribution function fitted to the velocity dispersion and 
line profile data. This is plotted as a function of the assumed potential's circular 
rotation velocity at the observed radius of the last kinematic data point. Filled 
symbols show best-estimate models derived with the optimal A = 10~ 4 ; open symbols 
represent models derived with A= 10~ 5 . The self-consistent (Af/L = const.) and the 
v c = const, models are marked separately. The horizontal dashed line shows the 95% 
confidence boundary derived from Fig.M 

These halo mass models have a constant density core, and the parameters are 
chosen such that the halo matters only at large radii. This is similar to the 
maximum disk hypothesis in spiral galaxies. We therefore call these luminous 
plus dark matter models maximum stellar mass models; we show a sequence of 
such models in Fig. ||. 

Dynamical models for NGC 6703 were derived from the data in a variety of 
these maximum stellar mass potentials, including those shown in Fig. [|. Because 
the dispersion profile is falling, only one of the two halo parameters can be 
determined. We have found that a well-constrained parameter is the circular 
orbit velocity at the radius of the last kinematic data point, v c (R m ). Fig. [| shows 
the results. A self-consistent potential is clearly inconsistent with the data, while 
a model with a completely flat rotation curve is marginally consistent with the 
data at ~ 1.5er. All models that are not ruled out at 95% confidence (according 
to Fig. |) have v c (R m ) = 250±40km S - 1 at R m = 2.6R e . This corresponds to a 
total mass in NGC 6703 inside 78" (13.5 kpc, where h 50 = # /50km/s/Mpc) 
of 1.6 — 2.6 x lO 11 /^ 1 M.Q, and to an integrated B-band mass-to-light ratio out 
to this radius of Tb = 5.3 — 10, rising from the central = 3.3 by at least a 
factor of 1.6. 

Fig. |6] shows the predicted kinematics from several models together with the 
data, and the implied anisotropy profiles. The anisotropy of the stellar distri- 
bution function in NGC 6703 thus changes from near-isotropic at the centre 
to slightly radially anisotropic (/? = 0.3 - 0.4 at 30", (3 = 0.2- 0.4 at 60"). It 
is not well-constrained at the outer edge of the data, where (5 = —0.5 — h-0.4, 
depending on variations of the potential in the allowed range. Notice that large 
anisotropy gradients between 40" and 70" and across the outer data boundary 
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Figure 6. Dynamical models for the kinematics of NGC 6703 in several luminous 
plus dark matter potentials, compared to projected velocity dispersion (top left 
panel) and VP-shape parameter h A (bottom left). The right panel shows the models' 
intrinsic anisotropy parameter [3(r), with the same linestyles: self-consistent model 
(stars only; long-dashed), v c = const, model (short-dashed), three models with 
v c (78") in the lower part of the acceptable range, according to Fig. (full), and two 
models with v c (78") in the upper part of this range (dotted lines). 

are implied if the rotation curve of NGC 6703 were completely flat. The model 
with v c = const is thus not a very plausible one for this galaxy. 

Rix et al. (1997) have analyzed the velocity profiles of the EO galaxy NGC 2434 
with a linear orbit superposition method. This galaxy provides an interesting 
contrast to NGC 6703 because it has an essentially flat dispersion profile. Its 
kinematics are likewise inconsistent with a constant-M/L potential, but are well- 
fit by a model with v c = const. This can be interpreted as a maximum stellar 
mass model in the sense defined above, in which the luminous component with 
maximal T contributes most of the mass inside R e . In this case, T ~ 6.3. The 
kinematics of NGC 2434 are also well-fit by a range of specific, cosmologically 
motivated mass models (Navarro, Frenk & White 1997) which, if applicable, 
would imply lower T, down to 4.3, and significant dark mass inside R e . 

It is clear that the dynamical analysis determines only the distribution of total 
mass, but cannot fix the mass-to-light ratio of the stellar component separately. 
To do this, we must, like in spiral galaxies, use extra constraints such as from 
population analysis. The fact that the central mass-to-light ratio in NGC 6703 
comes out low (3.3) even in a maximum stellar mass model may mean that the 
contribution of the halo to the mass within 0.5 — lR e cannot be very large. This 
needs to be studied in more detail; if true it would give an interesting constraint 
on galaxy formation theories. 
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5. Conclusions and Discussion 



The principal result of the work described here is that both the anisotropy 
and the mass distribution of an elliptical galaxy out to 2 — 3R e can be derived 
from velocity dispersion and line profile shape measurements. The two ellipticals 
analyzed so far do have dark halos, the resulting circular rotation velocity curves 
are nearly flat, and the mass-to-light ratios of the stellar populations have come 
out low. In NGC 6703 there are indications for deviations from a flat rotation 
curve in the sense expected if the luminous matter has segregated to the center; 
these need to be confirmed with more accurate observations. Future studies 
must show whether there are differences between the shapes of the true circular 
velocity curves of elliptical galaxies. 

The strength of the dynamical analysis of elliptical halos is that it simulta- 
neously yields the stellar-dynamical anisotropy structure, a further strong con- 
straint on formation theories. Moreover, since the dynamical analysis works best 
at high surface brightness, it is best-suited for constraining small halo core radii. 
The drawback of the dynamical analysis is that it can hardly be pushed further 
than to about 3R e , and that there are problems in deprojecting non-spherical 
systems (Gerhard and Binney 1996). These can, however, be eased by concen- 
trating on the flattest elliptical galaxies. For the E0 galaxy NGC 6703, current 
data determine the true circular velocity at the radius of the last kinematic data 
point to about ±40kms _1 at 95% confidence, within the framework of spherical 
models. The largest part of the quoted uncertainty is due to the finite radial 
extent of the data. 

X-ray observations have the advantage that they can yield the mass distribu- 
tion out to very large radii. Only few systems have been studied sofar, however, 
some of which are cluster brightest galaxies. A typical resolution limit is ~ R e . 
Gravitational lensing observations can either give an accurate determination of 
the mass within an Einstein radius (in the case of strong lensing), or a statisti- 
cal determination of the outer mass profiles of distant ellipticals (in the case of 
weak lensing). Lensing observations only test projected masses, but they have 
the advantage that dynamical equilibrium need not be assumed. 

We will therefore learn most about ellipticals by combining the different tech- 
niques, emphasizing their strengths and minimizing the uncertainties, and cover- 
ing as large a radial range as possible. For example, if the gravitational potential 
of an elliptical galaxy at several R e is known from X-ray data, this will sub- 
stantially reduce the error in the dynamically derived outer anisotropy. Stellar 
dynamical models derived from absorption line kinematics can be combined with 
discrete velocity samples at larger radii, such as from planetary nebulae searches 
(Arnaboldi & Freeman 1997), and stellar dynamical mass-to-light ratios are re- 
quired for normalisation in gravitational lensing studies (Kochanek & Keeton 
1997). 

The fact that also ellipticals have dark halos has been clear from X-ray and 
other observations and is not too surprising. However, the detailed study of their 
dark matter halos is of particular interest, since at least cluster ellipticals are 
now believed to be very old (see Bender et al. 1998 and references therein) and 
must have formed when the mean density of the Universe was still high. It will 
also be interesting to see whether the uniformity of ellipticals in the fundamental 
plane extends to their detailed dynamical properties, and whether this is caused 
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by dynamical processes (Merritt 1997) or reflects the details of the hierarchical 
formation process. 
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